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Nomenclature 

Uc global heat loss coefficient (W/m? K) 

Ug bottom heat loss coefficient (W/m? K) 

Ug edge heat loss coefficient (W/m? K) 

Ur top heat loss coefficient (W/m? K) 

hri i+1 radiation coefficient between two adjacent layers 
(W/m? K) 

hCii+1 convection coefficient between two adjacent layers 
(W/m? K) 

he radiation coefficient between absorber plate and 


first cover (W/m? K) 
V wind speed (m/s) 
I solar radiation (W/m?) 
Ta ambient temperature (K) 
Tp absorber plate temperature (K) 
Th water inlet temperature (K) 


Tio water outlet temperature (K) 

Te cover temperature (K) 

Toi temperature of the ith cover (K) 

Thri average temperature of the enclosed space (K) 

Ts sky temperature (K) 

E; emissivity of the ith layer (dimensionless) 

oO Stefan—-Boltzman constant (W/m? K~*) 

N number of glass covers (dimensionless) 

Nu Nusselt number (dimensionless) 

k conductivity of the air gap (W/m K) 

l plate cover spacing (m) 

Ra Rayleigh number (dimensionless) 

Pr Prandtl number (dimensionless) 

v kinematic viscosity of air (m° s7!) 

Ta optical efficiency of the collector (dimensionless) 

te transmission coefficient of the cover (dimension- 
less) 

Ap absorptance of the absorber plate (dimensionless) 

Dc reflection coefficient (dimensionless) 

hyn+1 wind convection coefficient (W/m? K) 

F collector efficiency factor (dimensionless) 

Cor specific heat of water (J/kg K) 

MM water mass flow rate (kg/s) 

n collector efficiency (dimensionless) 

D tube diameter (m) 

w width of the absorber area (m) 

B tilt angle (°) 

Ac collector area (m?) 


1. Introduction 


Nowadays, renewable energy is considered as the key to a 
sustainable energy future. It can have a beneficial impact on the 
environmental, economic, and political issues of the world. Fossil 
fuel price increase, climate change, adverse environmental impacts 
- in brief the respect of the Kyoto agreement - makes the 
exploration for a sustainable way to use energy more important 
than ever. Soteris [1] has recently shown that the benefits arising 
from the installation and operation of renewable energy systems 
can be classified into three categories: energy saving, generation of 
new working posts, and the decrease of environmental pollution. 
Solar energy may play an outstanding role to substitute 


non-renewable energy sources as it may also play an important 
prerequisite to enhance sustainable development. 

Solar water heating systems (SWHS) are very common systems, 
extensively used in many countries with good solar radiation 
potential, such as in the Mediterranean countries. SWHS are often 
viable to replace electricity and fossil fuels used for many home 
applications [2]. In general in a SWHS, conventional flat plate solar 
collectors (FPSCs) with a metal absorber plate and covers are used 
to transform solar energy into heat [3]. In a FPSC, the incident solar 
radiation is converted into heat and transmitted to a transfer 
medium such as water. 

Several FPSC designs have been proposed and discussed in 
literature. Different aspects can be investigated comparing the 
performances of those systems. These aspects can be mainly 
classified in: the thermal performance, lifetime and cost, 
sustainability or durability, maintenance and installation devices 
[4]. 

FPSC thermal performances can be modelled according to 
complex functions of the meteorological data, of the design 
system and finally of the operating conditions. Since the 
meteorological conditions cannot be modified, the FPSC per- 
formance can be enhanced working either on the operations 
conditions or on the SWHS components. In the design of SWHS 
components, care should be especially given to the number and 
nature of the transparent covers and on the selectivity of the 
absorber [5]. Glass is a quite common choice as a cover for solar 
thermal devices since it absorbs almost all the infrared radiation 
(IR) re-emitted by the absorber plate, and then leads to a 
reduction of the greenhouse effect which is responsible of the 
enhancement of the thermal efficiency of the solar collector. On 
the other hand, the use of a glass cover in rural zones of 
developing countries has two major disadvantages: its high 
installation cost, and its fragility both during transportation and 
in service [6]. So, some authors [7,8] recommend the use of 
plastic covers because of their relatively low weight and good 
resistance against shocks. On the other hand, the long-term 
plastic cover exposition leads to degradation. 

Although theoretical simplified approaches have been avail- 
able for many decades, an exhaustive FPSC thermal analysis is still 
a hard task since several characterizing parameters are required. 
In literature, there are numerous studies on the SWHS perform- 
ance analysis. 

Akhtar and Mullick [9] proposed a model to compute the FPSC 
glass cover temperature with single glazing, and then [10] 
enhanced this model proposing a set of relations to compute 
the glass cover temperatures for double glazing. The values of top 
heat loss coefficient are computed and are close to the numerical 
solutions of heat balance equations. 

Smyth et al. [11] discussed integrated collector storage solar 
water heaters. Kaldellis et al. [12] presented an analysis of the time 
variation of the local market taking seriously into account the 
domestic solar water heating system annual replacement rate. 
Batidzirai et al. [13] demonstrated the importance and significance 
of water heating in the national energy end-use regime in 
Zimbabwe and concluded that solar water heating can result in 
economic, social and environmental benefits for the whole 
country. It shows that water heating in the residential, health 
and tourism sectors contributes about 13% of the coincident winter 
peak demand and accounts for 8% of final electricity consumed in 
the country. Muneer et al. [14] investigate the possibility of 
exploitation of solar energy within the Turkish textile industry 
with particular reference to thermal application and presented a 
detailed life cycle assessment and relevant economics of solar 
water heater. Muneer et al. [15] proposed an introduction of built- 
in-storage water heaters for Pakistani textile industry. They 
presented a comparison of the thermal performance of two 
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different designs of built-in storage water heater plain and newly 
designed finned type. Armando and Oliveira [16] presented a new 
approach to express the long-term performance of general solar 
thermal systems. Jannot and Coulibaly [6] presented a new set of 
equations for the radiative balances of the absorber plate and the 
transparent cover of a solar air heater covered with a plastic film. 
Janjai et al. [3] developed a mathematical model to simulate the 
performance of a large area plastic solar collector. Bilgen and 
Bakeka [17] designed and studied a simple solar system and 
evaluated its economics and thermal performance. Soteris [1] 
studied the thermal performance, economics and environmental 
protection offered by thermosiphon solar water heating system 
and showed that a considerable percentage of the hot water needs 
of the family can be satisfied by solar energy. 

The general aim of this work is to support decision makers in 
the FPSC design phase and in the control of the inlet water flow 
in the FPSC in order to optimize the FPSC efficiency. This aim is 
achieved according to a two-folded approach. Firstly, as 
described in Section 2, a model is formulated and implemented 
taking into account various modes of heat transfer in the 
collector. The goal is to investigate the impact of the number 
and of the types of the covers on the top heat loss and the 
related thermal performance in order to support decision 
makers about the most cost-effective design. The proposed 
model can also be used to investigate the effect of the different 
parameters which may affect the performance of the collector. 
Secondly, as described in Section 3, a two objective constrained 
optimization problem has been formulated and implemented to 
evaluate the optimality of the performance of different design 
approaches. The goal is to support decision makers in the 
definition of design parameters, i.e., the absorber plate area, the 
pipes’ diameter, the mass flow rate, and distance between the 
tubes. Section 4 shows an application of the proposed 
methodology with environmental temperature and irradiation 
characteristics of Tétouan (Morocco). Conclusions and future 
research directions are finally shown. 


2. FPSC heat transfer modeling 


The evaluation of the thermal loss coefficients is the funda- 
mental task to assess the FPSC performance. 

The FPSC global heat loss coefficient Ug (W/m? K) is the sum of 
the top loss (Ur), the bottom (Ug), and the edges (Uç) coefficients. 
This loss value is established between the collector and its 
surrounding by conduction, infrared solar radiation and convec- 
tion [18]: 


Uc = Ur + Ug + UE (1) 


According to Ref. [19], in an efficient design, the edge and the 
bottom loss coefficients are smaller than the top heat loss 
coefficient. In addition, assuming that the FPSC insulation of the 
edges and of the bottom is responsible for the resistance to heat 
flow, the sum of the edge and bottom heat loss coefficients can be 
set to a constant value A. 


Uc = Ur +A (2) 


In this paper, the top heat loss coefficient is referred to a FPSC 
model with N covers. Specifically, the heat loss is referred to 
several FPSC layers, identified by the index i, i = 0,. ..,N + 1, where 
the index 0 is referred to the plate adsorber p, i= 1 is referred to 
the bottom cover, i = 2 is referred to the second cover (if present), 
i=N is the top cover, and N+ 1 is the external environment or 
ambient a. 


The overall energy loss coefficient through the top, Ur, can be 
assumed as the result of convection and radiation heat losses 
between parallel plates (or layers). That is: 


N-1 a 
Uy = (Sere + mran) ') i=0,...,N-1 (3) 


i=l 
The radiation coefficients hrjj+1 (W/m? K) are expressed by Ref. [3]: 


o(T; + Tix1)(T? + T?.,) en ee (4) 


Arig = (1/Ei) + (1/Ein1) - 1 


oEn(Ty — TS) 


Ty T; (Abis) 


hrn N41 = 
where hry,n+1 represents the radiation exchange between the 
outer cover and the sky, Ey (dimensionless) is the emissivity of the 
Nth cover, Ty (K) is the temperature of the Nth cover (cover with 
direct contact with the ambiance), T, represents the ambient 
temperature (K) and T, is the considered to be the sky temperature 
(K). 

The convection heat loss coefficients hcii (W/m? K) are 
conducted through the free convection heat transfer existing in 
the incarcerate air between the adjacent surfaces. Eq. (5) defines 
the hcii coefficient between the absorber and the first inner 
cover, as well as between two adjacent covers [5], as a function of 
kii+ı (W/m K), the conductivity of the air between the layer i and 
i+1, and the distance lii (m) between the layer i and i+ 1. 
Eq. (5bis) empirically defines the outer convection loss coefficient 
in term of the wind velocity V (m/s). 


Ki isa Nui i41 


hCi i41 = =0,..., N-1 (5) 


liir 


hy N41 =2.8 +3V (5bis) 


Regarding the experimental investigations in the literature [20], 
the relationship between the Nusselt number and Rayleigh 
number for tilt angles of the FPSC from 0° to 75° is: 


Nui = 1 
i 16 x 
1.44|1 1708(sin 1.88) 1 1708 
Rajji., cos B Raji, cos B 
+ 
Rajj,1 cos A\ 
( a ‘| i=0,...,N—1 (6) 


where the apex (+) is referred to the condition that only positive 
values of the terms in the square brackets are used. 
The Rayleigh number Raq; (dimensionless) is given by Ref. [20]: 
Prii (Ti — Tia), 
8Priiz1 (Ti i1) Wt Gig N—1 (7) 
Tm.v? 


m-"iit+1 


Rai i1 = 


where for i= 0, To = Tp, and Prii+ı (dimensionless) is the Prandtl 
number at the air gap between the layer i and i+1 which 
depends on the average temperature Tmii+ı = (Ti+ Ti + 1)/2; g 
(m s7?) the gravitational acceleration, and v;i (m° s7!) is the 
cinematic viscosity at the air gap between the layer i and i+ 1. 


2.1. Computational approach 


The top heat loss may be computed according to several 
empirical relations. Hereinafter, an algorithm is proposed to 
compute it, on the basis of the direct computation of each cover 
temperature. Considering a steady state energy transfer, the energy 
exchange from one gap air to the other is equal to the overall energy 
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exchange between the absorber and the ambient. Under this 
hypothesis, the following system of equations can be defined: 


Ur(Tp — Ta) 
(hep + hrp1) 
id (Arp + hcp1)(Tp = T1) 


T =Tp 


T2 


hci 2 + hr12 
T-T (Arj_2i-1 + NCi_2 4-1) (Ti-2 — Ti-1) 
aes heat hrii (8) 
Ty =T] 4 (hr y + he)(Tn-2 — Ty-1) 


hcy-1,N + hrN-1,N 


In the case of one cover, taking into account the first equation of the 
system of Eqs. (8), the equation defined by Eq. (3), and two 
unknown variables Ur and T;, a solution can be found in an 
iterative two-step approach as described below. 

In the first step of the first iteration (i=1), the ambient 
temperature T, and the temperature of the plate T, are known. 
In addition, an hypothesis on the temperature of the first cover, 
T!, is made, where the index 1 is referred to the cover and the 
apex 1 is referred to the iteration. Specifically, T! is supposed to 
vary in a range T! € [Ta,Tp], for example taking into account an 
average value T} = (Ta + Tp)/2. From Eq. (3), it is possible to 
compute Ur. 

In the second step, the current value of Ur, as well as the T, and 
Tp values, are given as input to the first equation of (8), anda new T; 
is obtained. 

A convergence test is so evaluated. If Ir E f | <é, where £ is a 
suitable small value such as ¢=10°°, the algorithm ends. 
Otherwise a new iteration i=2 is made. In this new iteration, 
Ti = Ñ. This process will be iterated until the convergence test is 
satisfied. 

If other covers are presents, their temperature can be computed 
using the same process used for the case of one cover. 


2.2. Reference models 


In the case study, the results given by the model described 
above are compared with the results given by two different 
methods, Chaudhuri [21] and Akhtar and Mullick [9,10]. These 
methods allow the computation of the top FPSC heat loss 
coefficient. A description of these methods is briefly quoted 
hereinafter. 


2.2.1. Klein method 
The top loss coefficient Ur can be computed using a conventional 
empirical equation proposed by Ref. [21], and which is: 


Ta N 11` 
T |(849/Tp) [Tp —Ta/N+ 9002 finns 
o(T3 +12) (Tp — Ta) 
(Ep + 0.05N(1 — Ep)! + [(2N + — 1)/Eg] — N 


(9) 


where $= (1 — 0.04hn n+ + 0.005h%,y,,,)(1+0.091N), N is the 
number of glass covers, E, (dimensionless) is the emittance of the 
plate, and E, (dimensionless) is the emittance of glass. 


2.2.2. Akhtar method 

In the Klein empirical relation, the global top heat loss is 
evaluated since the cover temperature is unidentified. Thus, the 
condensed regrouping parameters in the Klein relationship are 


deficient in the determination of the radiation and convection heat 
losses coefficients. For these reasons, and to make our approach 
more consistent, the approximation method of Akhtar and Mullick 
[9,10] is more appropriate to compare the obtained results. 
Actually, this method aims to provide an empirical relation for the 
glass cover temperature T, while bearing in mind the heat transfer 
processes in the FPSC. A more detailed description of the method 
can be found in Akhtar and Mullick [9,10]. 


3. Optimization of the FPSC design 
3.1. FPSC thermal performance 


The solar energy, intercepted by the plate absorber is equal to 
the incident solar radiation reduced by optical losses, which is a 
function of the optical efficiency. The FPSC optical efficiency 
strongly depends on the used material type, such as transparent 
covers. The cover characteristics influence the energy transferred 
to the absorber, and consequently the thermal performance of the 
overall system. So, the selection of a suitable component constitutes 
a fundamental task. The optical efficiency or the effective product 
transmittance-absorptance is expressed by Ref. [5]: 


_ TcQp 
STET 
where (tæ) is the FPSC optical efficiency (dimensionless), Te is the 
cover transmission coefficient (dimensionless), œp is the absorber 
plate absorptance (dimensionless), and p: is the reflection 
coefficient (dimensionless).The FPSC performance is described 
by an energy balance that converts the incident solar energy into 
useful energy gain, thermal losses and optical losses. Assuming a 
steady state, the energy balance of different components is 
expressed by: 


(10) 


. . dT j 
MieCpp -gy = FUTAI — Uc (Tr — Ta)] (11) 
with 
r 1/Uc (12) 


~ w[(1/Uc[D + (w — D)F]) + (1/7Dhr;)] 


where mr (kg/s) is the water flow in the FPSC tube, T;, is the outlet 
water temperature (K), cpf (J/kg K) is the specific heat of water, F 
(dimensionless) expresses the FPSC efficiency factor, D (m) is the 
water tube diameter, w (m) is the width between two centres of 
adjacent tubes, hg (W/m? K) is the radiation coefficient between 
absorber plate and first cover and F is the standard fin efficiency 
(dimensionless) that depends on D, w and mg. 

Assuming that F and Ug are independent from the x position 
(that is referred to the length of the water tube), the solution of 
(11) is: 


Tio — Ta- (ta)I/Uc 
Ti — Ta — (ta)I/Ug 


= exp(—ULF’Ac/MC pF) (13) 


where Ta is the inlet water temperature and A, (m7) is the area of 
the FPSC. 

In the computations quoted hereinafter, the sum of the bottom 
and the edges losses coefficient A is supposed to be constant and 
equal to 1.2W/m?K. A measure of the instantaneous FPSC 
performance is the collector effectiveness, expressed as the ratio 
of the instantaneous useful thermal power extracted from the 
water collector Q,(t) (J) to the instantaneous incident solar 
irradiation A.l(t) (J), and it is defined as reported by Farahat et al. 
[22]: 


n(t) ere (14) 
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The evaluation of the daily FPSC efficiency takes place by 
integrating the previous equations over all the operating solar 
system hours. This expression could be obtained by integrating 
from the instant t, of sunrise to the instant t, of sunset. 

The overall thermal efficiency ņ during the daytime is 
calculated from the relationship: 


7 ie MeCpe (Tho — Tr )dt 
Ac f I(t)dt 


(15) 


where 


ts 


MpCpf(Tf> — TR )dt 


Qu(t) = (15bis) 


tr 


3.2. The optimization problem 


The model described in Section 2 highlights the importance of 
choosing proper design parameters, such as the sizes of the 
absorber plate, the pipes’ diameter, etc., for the performance of the 
heating through a FPSC. In this section, an optimization approach is 
followed in order to ensure an optimal performance of the FPSC 
design. Specifically, the optimization model defined hereinafter 
aims to find the optimal design of the solar system that offers a 
maximum efficiency and a maximum outlet water temperature. 
So, the problem is formalized as a two-objective decision problem, 
and its application may be useful for decision makers, as it gives 
the opportunity to obtain the optimal value of several important 
design options. Moreover, the Pareto optimal solutions, as shown 
in the case study, may be useful for decision makers since it is 
possible to identify all the possible scenarios, and, depending on 
the available financial resources, either a specific solar collector 
area or a particular working fluid mass flow can be chosen. 

In the proposed optimization problem, the decisional variables 
are related to the main design parameters, i.e., the absorber plate 
area A. (m7), the pipes’ diameter D (m), the mass flow rate m; (kg/ 
s), and distance between the tubes of the FPSC, w (m). Two non- 
commensurable objectives have to be maximized: the outlet 
temperature of the water, Tto (K), and the efficiency of the overall 
collector, 7 (dimensionless). 

In a discretized time, the overall objective function is given by 
the weighted sum of the two contributions. That is: 


ts 
Z= SCT +N (16) 
t=tr 
where g (K) is a weighting factor and the optimization instants t 
are given in hours. 
Based on the continuous formulation of Eq. (13), the outlet 
temperature of the water is given by: 


(tæ) Ta — T! — (tæ) 


Tt, (Ac, w, D, rip) [r ' 


Ug ` Ug 
7 
x exp UEA =1,...,T (17) 
MfCpf 


where the apex tis the time instant, at which ambient temperature 
and irradiance measurements are also available. 

Based on the continuous formulation in Eq. (15), the overall 
efficiency is: 


o =I MeCpe (Th, (Ac, w, D, mir) — Ti) 
Ac ye É 


rt=1 


n (18) 

where t= 1 is the time of sunrise and t=T is the sunset time. 
Different constraints have been taken into account in the 

optimization problem and are related to the size that may assume 
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Fig. 1. Climate data used in the simulations. 


the different decision variables. Specifically, there are constraints 
that limit the size of the diameter, the area, and the distance among 
tubes to a range of minimum and maximum values. That is: 


Dmin = D < Dmax 
Wmin < W < Wmax (19) 
Amin < Ac < Amax 


where Dyin, Dmax, Wmin, Wmax, Amin, Amax are known parameters. 

Finally, there are constraints related to the operating param- 
eters that require the existence of maximum and minimum values 
for the flow rate of the fluid mass. That is: 


Mmin < my < Mmax 
where Mmin and Mmax are known parameters. 
4. Results and discussions 


Climatic data that are used in the current investigation 
correspond to the case of a clear type day of July, measured in 
Tetouan (35°35’N 5°23’W), Morocco [23]. The hourly diffuse 
radiation has been calculated from the hourly global radiation, 
using a relation available in the literature [23]. The climatic data 
are reported in Fig. 1. 
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Fig. 2. Top heat loss coefficient versus the absorber plate temperature for one 
Plexiglas cover. 
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Fig. 3. Top heat loss coefficient versus the absorber temperature for double Plexiglas 
covers. 


4.1. Modelling FPSC heat losses 


The top heat loss coefficient Ur is calculated for the FPSC using 
the iterative numerical method, which seems to be among the 
most accurate one. A value of 107° was considered as a residual 
between two consecutive values of cover temperatures. In order to 
assess the validity of the proposed method to predict the top heat 
loss coefficient, the results are compared with some of those 
existing in the literature. Simulations were first performed for a 
single cover, at given values of the absorber emissivity, 0.1 and 
0.95, which are considered respectively as minimum and maxi- 
mum values. In Fig. 2,a good agreement was observed between our 
results and those given by Akhtar and Klein methods. Figs. 3 and 4 
display the top heat loss coefficient versus the absorber tempera- 
ture, for double covers; thus respectively for a double Plexiglas 
covers, and combined glass Plexiglas covers. 

In order to evaluate the performance of the FPSC, Fig. 5 displays 
the top heat loss coefficient as a function of the absorber 
temperature, for the three test cases. The heat losses are 
diminished for a selective plate absorber (Ep =0.1), this kind of 
cover re-emits a small part of the solar incident irradiation. It can 
be also remarked that for a selective absorber, either for double 
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Fig. 4. Top heat loss coefficient versus the absorber temperature for combining 
glass—Plexiglas covers. 
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Fig. 5. Top heat loss coefficient versus the plaque absorber temperature for one, and 
two covers (glass/Plexiglas and Plexiglas/Plexiglas). 


Plexiglas covers or combined glass Plexiglas covers, the top heat 
loss coefficient is approximately equal. However, for a non- 
selective absorber, double plastic covers were those that give us a 
slightly higher heat top losses coefficients. 

Hence, the combined plastic—glass covers were thus judged to 
be the most suitable ones in terms of reducing the top heat loss 
coefficient. Fig. 6 displays the comparison between the ratio of the 
top heat loss coefficient and the top heat loss at 45°, as a function of 
a given tilt angle of the collector system. In this latter illustration, 
results are presented in terms of a normalized value of the heat loss 
similar to those given by Ref. [19]. It appears that the current 
results are similar to those presented by Akhtar and Mullick [9,10], 
and this for both values of emissivities. 

The FPSC performance depends on different design compo- 
nents. This study focuses on an assessment of the effect of the 
absorber emissivity on the heat transfer process. Fig. 7 displays the 
heat loss coefficient as a function of the emissivity, and that for 
different type and number of covers. It can be recognized that, 
whatever the emissivity value, a single Plexiglas cover was the one 
that gives a higher value of Ur. Similar results are found for double 
covers (glass/Plexiglas and Plexiglas/Plexiglas). Similar trends are 
observed between the two kind of double covers. 
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Fig. 6. Ratio of the top heat loss coefficient and the top heat loss coefficient 
corresponding to 45° as a function of the tilt angle. 
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Fig. 7. Top heat loss coefficient against the emissivity of the absorber for single and 
double covers. 


As a second part, investigations were done on the effects of the 
variation of the FPSC operation conditions. These investigations 
included the effects of changing the operation inlet water tempera- 
ture ata fixed value of flow rate, and, then, the rate at which the water 
flows in the tubes at a given value of the inlet flow temperature. 

Studies were carried out by using two kind of double covers: the 
double plastic and the combined glass plastic covers. As displayed 
in Fig. 8, practically, there is no significant difference between the 
maximum outlet temperature reached by the flat plate for the 
three values of the inlet temperature, and thus Fig. 9 shows the 
same trend for the combined covers. 

In order to give more details on the performance of both types of 
covers, Fig. 10 shows the outlet collector temperature at a given 
value of water flow rate (mı f = 0.001 kg/s) and at given value of 
inlet fluid temperature (Ts = 20°). It appears that the combined 
covers provide a better outlet performance than the double plastic 
covers. A difference of 7° is observed. 

In general, the mass flow rate of the working fluid is one of the 
most important parameters affecting the performance of solar 
collectors. In this study, the influence of the mass flow rate on the 
water outlet temperature and the efficiency was studied. Fig. 11 
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Fig. 8. The effect of inlet water temperature on the hourly variations of the outlet 
temperature for a double Plexiglas cover. 
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Fig. 9. The effect of inlet water temperature on the hourly variations of the outlet 
temperature for a combined Plexiglas glass cover. 
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Fig. 10. Comparison of the hourly outlet water temperature for one, combined glass 
Plexiglas cover and double Plexiglas cover for Tg = 20° and m = 0.001 kg/s. 
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Fig. 11. Efficiency of the FPSC for different kinds of covers. 
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Fig. 12. Maximum outlet water temperature and efficiency of the FPSC versus the 
water flow rate. 
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Fig. 13. The collector efficiency versus the output average temperature of water. 


shows the effect of the water flow rate on the FPSC daily efficiency, 
for a single Plexiglas cover, double Plexiglas covers and combined 
glass Plexiglas covers. This comparison reveals that the glass- 
Plexiglas covers have much better behavior than the others one. 
This is especially due to the fact of combining the performances of 
glass, which has a better transparency to the incident irradiation 
and the Plexiglas devices that are good for collecting the 
greenhouse effect. The results in Fig. 12 show that the maximum 
peak of outlet temperature increases with the decrease of flow 
rate, while an opposite trend is observed for the effect on the 
collector efficiency. A value of 2.66 x 107° kg/s is considered to be 
suitable flow rate to maintain good efficiency and adequate fluid 
outlet temperature. 


4.2. The optimization results 


The necessary parameters to solve the optimization problem 
are reported in Table 1. In particular, the non-zero values of solar 
irradiation and outside temperature are reported. Instead, in Table 
2, the non-time dependent parameters are reported. 

The optimization problem is non-linear with continuous 
decision variables. It has been solved for a time horizon of 24 h 
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Fig. 14. The variation of collector efficiency and the daily average temperature of 
water as function of water flow rate. 
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Fig. 15. The variation of collector efficiency and the daily average temperature of 
water as function of collector area. 


using the optimization package Lingo 9.0 (www.lindosystems.org). 
It takes less than 1 min to obtain the optimal local solution.! 

Fig. 13 shows the results obtained from the optimization 
process as regards the behavior of the two terms of the objective 
function, i.e., the efficiency versus the mean outlet temperature. 
The evolution of the solutions as regards the water flow rate and 
the collector area according to the objective function are reported 
respectively in Figs. 14 and 15. Fig. 14 shows the variation of the 
average outlet water temperature and the collector efficiency 
according to different values of the weighting factor, thus fixing the 
value of the mass flow rate (at 0.0037, 0.0058, 0.0074, 0.0087 and 
0.01 kg/s). According to this figure, the trend of curve goes by 
convergence, in others words, the gap between the two set of 
temperature and the efficiency diminishes as the value of the mass 


1 The multistart optimisation option has been chosen to enhance the reliability of 
the local optimal solution that have been obtained. 
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Table 1 
Irradiance and ambient temperature of a clear type day of July in Tetouan. 
Time (h) 6 7 8 9 10 11 12 13 14 15 16 17 18 19 
T! 23.9 24.4 25 25.5 26 26.5 27 27.5 27.75 27.9 27.8 27.5 27.3 26.6 
I 0 150 375 600 750 850 875 875 850 750 600 400 175 0 
Table 2 oo glass covers gives a better performance to the FPSC. Results of the 
Inputs data for optimization. developed model have been compared with other ones available in 
Dmin=0.01 m Dmax= 0.05 m the literature, and similar results have been found. Finally, a two- 
ee m a objective non-linear optimization problem has been formalized and 
M =0.001 m3/s M,,...=0.008 m?/s solved though mathematical programming technique. Moreover, 
hg =300W/m?K m=5 the optimization problem has been run for different values of the 
(ta)=0.614 T= 20° weighting factor (that represents the relative importance given to 
= = 2 e : a . . F 
cpr=4190 Uc =5.31695 W/m" K the two objectives) in order to build the Pareto optimal solutions. 
The aim of this approach is to find different configurations of the 
100 60 system in order to help decisions makers to choose among a variety 
of design and system operations. The results show that the water 
f desig d syst perat Th Its show that th t 
ae 50 mass flow rate significantly influences the values of the two terms of 
S ð the objective function, whereas, the effect of the collector area does 
Š > not have a considerable effect on finding a better compromise 
3 70 40 X between the two conflicting objectives. 
8 = 
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Fig. 16. The effect of the weight factor on collector efficiency and daily average 
temperature of water. 


flow rate increases. Whereas, in Fig. 15, which shows the variation 
of the average outlet water temperature and the collector 
efficiency according to different values of the weighting factor, 
and fixing the values of the collector area at (4, 6, 8 and 10 m2). In 
this figure, a divergence of the gap between the average outlet 
water temperature and the collector efficiency is remarked. The 
results shown in these two figures are particularly helpful for the 
selection of the optimal FPSC design. The trends detect the 
behavior of the collector efficiency and they are associated to the 
daily average temperature and different decisional variables. Thus, 
the decision makers, according to their funds and the collector 
performance requirements, will select the suitable collector area 
and the operating water flow rate. They will envisage all the 
possible scenarios, and also choose among a variety of system 
designs. Fig. 16 shows the variation of the optimal objective 
function when varying the weighting factor. It can be seen that the 
efficiency increases with the weighting factor until a critical value 
of œ =2198 K, after which the efficiency does not vary and it is 
equal to 52.11%. The mean temperature decreases until the same 
critical value of 2198, and then remains constant at 56.42 °C. 


5. Conclusion 


In this paper, investigations of the thermal behavior of the FPSC 
have been done. The goal is to compare the overall loss coefficient of 
the collector under different numbers and kinds of class covers. The 
results depicts that the top heat loss (Ur) is lower when adding 
another cover to the FPSC, it was found that combining plexigals and 
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